Flux expulsion and greedy bosons: frustrated magnets at large 
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We investigate the Sp(Af) mean-field theory for frustrated quantum magnets. First, we establish 
some general properties of its solutions; in particular, for small spin we propose simple rules for 
determining the saddle points of optimal energy. We then apply these insights to the pyrochlore 
lattice. For spins on a single tetrahedron, we demonstrate a continuous ground state degeneracy 
for any value of the spin length. For the full pyrochlore lattice, this degeneracy translates to a 
large number of near-degenerate potential saddle points. Remarkably, it is impossible to construct 
a saddle point with the full symmetry of the Hamiltonian — at large A'^, the pyrochlore magnet 
cannot be a spin liquid. Nonetheless, for realistic finite values of N, tunnelling between the nearly 
degenerate saddle points could restore the full symmetry of the Hamiltonian. 



Introduction. The behavior of quantum magnets 
with strong frustrationi is one of the central open ques- 
tions in the study of modern magnetism. The most cele- 
brated members of this class of problems are the nearest- 
neighbor Heisenberg antiferromagnets on the kagome and 
pyrochlore lattices and the challenge is to work out their 
phases at varying values of the spin S and temperature 
T. At large S, semiclassical computations become fea- 
sible if not straightforward!^ Small spins have typically 
been treated by entirely separate methods, and for the 
case of the pyrochlore lattice these have involved starting 
from a symmetry-breaking decomposition of the lattice.'^ 

As a complement to this work, we employ a large- 
N method that can be applied over the full range of 
S and, particularly for the work described in this pa- 
per, allows us to keep all symmetries in the formalism 
at small S. Specifically, we consider the enlargement of 
the SU(2) = Sp(l) symmetry group of the Heisenberg 
model to Sp(A^). In the limit N —* oo, the problem re- 
duces to a mean-field theory for Schwinger bosonsr->^'5 
The mean-field theory can be improved by evaluating 
corrections. This method has been applied to the 
kagome lattice^, and found to predict a robust selection 
of a magnetically ordered state at large S melting into a 
disordered spin-liquid state at small S. 

We first establish some general properties of mean-field 
ground states of Sp{N) magnets on regular lattices at 
small values of spin. The main variables of the mean- 
field theory are complex numbers Qij defined on links 
{ij) of the lattice and representing probability ampli- 
tudes of finding a valence-bond singlet. We derive a per- 
turbative loop expansion, wherein spin length is used as 
a small parameter. To the lowest nontrivial order, va- 
lence bonds are found exclusively on links with largest 
exchange coupling. At higher orders, we find that a U(l) 
flux, constructed from phases of the valence-bond am- 
plitudes, tends to be expelled. These observations sys- 



tematize previous results obtained for several quantum 
antiferromagnets in the large- framework^ 

We then discuss the case of the pyrochlore antiferro- 
magnet where we report three principal results: (i) For 
a single tetrahedron, the two-parameter degeneracy of 
classical spins persists for any spin length despite the 
quantum effects encoded in the Sp(A^) computation! (ii) 
All of these solutions break a symmetry — spatial or time 
reversal — and exhibit bond or chiral order, (iii) Embed- 
ding such states in the pyrochlore lattice results in a 
vast number of solutions whose energies are very close to 
one another at small spin. This behavior is fundamen- 
tally different from the triangular and kagome lattices, 
where the principle of flux expulsion alone was sufficient 
to fix the ground states. The import of (ii) and (iii) is 
that, if the pyrochlore lattice is indeed a spin liquid at 
iV = 1, this will necessarily require tunnelling between 
saddle points. Such a scenario is not at all unlikely. 

Sp(A^) Hamiltonian. The Hamiltonian of an Sp(A^) 
antifcrromagnct may be written as^ 

H = -{1/2N) J2 J^j4j-^'^j, = e^0b,^b,p. (1) 

The antisymmetric tensor Sap is a block-diagonal 2Nx2N 
matrix 







(2) 



and (Ty is the 2x2 Pauli matrix. Here, bia is the anni- 
hilation operator of a boson of species a at site i. For 
N = 1, Eq. ((TJ reduces to the SU(2) Heisenberg Hamil- 
tonian written in terms of Schwinger bosons, related to 
spin operators via S'f = ^la'^ap^ipl'^ O'^^'^ doubly 

repeated fiavor indices is implied); the number of bosons 
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determines spin length: bj^bia = 2S. In the large- gen- A scaling transformation Qi 



X —^ a\ does not 



eralization, spin length is related to the number of bosons 
per flavoi^ k = bj^bia/N. 

Mean-field approximation The Sp(A^) mean-field 
equations,^ which become exact for N oo, involve a 
decoupling of the quartic terms in Eq. |^ with the aid 
of the link variables (6^6^73) = Q,j£a/3/2: 

,t 



i^MF = {J^J/2)J2{N\Q^J\^ - Q^Jec,pblJ^ 



Qtj^apbiabjp) + ^ Xi(b\j)^a " nN). (3) 



The chemical potential keeps the average number of 
bosons fixed at k,N on every site. The mean-field equa- 
tions are d{HMF) /dXi = (constraints on the boson 
numbers) and d{Hii-p) / dQij — (minimization of en- 
ergy). 

The mean-field Hamiltonian ((Jl is a sum of N identical 
copies, each containing two flavors only (up and down in 
the Schwinger-boson language). The energy of each copy 
is 0(1) in N . Therefore, two different vacua will have an 
energy difference 0{N), i.e., well separated in the limit 
N ~i 00 considered here. The low-lying excitations are 
5 = 1/2 bosons (spinons) whose energy is 0(1). 

For a total of TV sites, this Hamiltonian can, for each 
pair of flavors, be written in terms of TV-component row 
vectors and b^, column vectors bi^ and and TVxTV 
matrices Pij — JijQij/2 and Ajj = XiSij. A Bogoliubov 
transformation diagonalizes the part quadratic in bosons 
and yields a diagonal matrix of eigenfrequencies Q for the 
bosonic spinons giving the energy per flavor 



{Hmf)/N = Tr [PQV2 - (k + 1)A + n] 



(4) 



Uniform A. In the remainder of this paper we con- 
sider the mean-field theory with two restrictions. Firstly, 
we have tacitly assumed the absence of a condensate of 
bosons, (bia) = 0. (This translates into the lack of mag- 
netic order, (Si) — 0, for SU(2) spins.) This regime, 
dominated by quantum fluctuations, always exists for 
small values of spin length k.^'^ In addition, we will only 
consider states with uniform chemical potential, so that 
A = Al, where 1 is the TV x TV unit matrix. This simpli- 
fies calculations as matrices A and P commute. 

With said restrictions, the boson spectrum ri„ can be 
obtained from the eigenvalues v'^ of the matrix PP^: 



= - 1^1 det (PPt _ j,2i) = 0. 
The expectation value of the Hamiltonian is then 



(5) 



{Hmf)/N = Tr PQV2 + VA^I - PPt - A(l + k)1 
The boson-number constraint gives the condition 

l-HK = Trf(l-PPVA^)"^/^l /Trl. (6) 



affect the constraint equation Minimization of the 
energy determines the optimal scale a and yields 



N 



{Tr [(PPVA') (1 - PPt/A2)"'/'] } 



2Tr(POt/A2) 



(7) 



Further minimization of the vacuum energy IJJ) is done 
by varying the relative strengths and phases of the link 
variables Qij, subject to constraint l^. Our ad-hoc as- 
sumption of a uniform chemical potential restricts the 
choice of trial states in that all sites must be equivalent. 

Loop expansion at small k. The structure of 
the mean-field equations © and {Tj) suggests a solution 
by expansion in powers of PPt/A^. Taking the trace 
makes it a loop expansion: a generic term of the Taylor- 
expanded right-hand side of Eq. © has the form 



Tr(PPt) =^P„,(_p;j 



P 



(8) 



where we have accounted for antisymmetry, Pcb = —Pbc- 
Expressions for k and E = (H) involve sums over all pos- 
sible closed loops abc . . . za oi even length. The U(l) flux 
defined in Eq. ((SJ will play an important role below 
Convergence is particularly good when spinons have 
a short correlation length, as is the case for small k. In 
this limit, the physics is determined by the shortest loops. 
Formally, the loop expansion is a series in powers of k. 
In the following paragraphs, we develop the loop expan- 
sion and demonstrate that it leads to simple organizing 
principles for the behavior of Sp(A^) antiferromagnets in 
the quantum limit of small k. A variant of this strategy 
at high temperatures has been described previously^ 

Shortest loops, greedy bosons. The lowest order 
in PPt/A^ yields the energy per site, per flavor: 



{HI 

NN 



Tr (PPt) 
Tr(PQt) 



2 J2{ij) JijlQijl'^ 



(9) 



To leading order ("loops" of length 2), the energy de- 
pends on the absolute values, but not the phases of the 
link variables Qij. Minimization can be easily done if 
one interprets Jy P as a probability distribution. The 
energy is then simply the expectation value of — Kjy/2. 
An optimal probability distribution then will have zero 
probabilities Jij\Qij\'^ for all links except those with the 
largest Jy . For example, on a square lattice with first 
and second- neighbor couplings Ji,J2 > 0, the second- 
neighbor bonds will vanish, Q2 = 0, if Ji > J2; similarly, 
Qi = if Ji < J2 in this approximation. Hence 

Theorem (greedy bosons): in the hmit of 
small K, bosons form bonds Qij 7^ on the 
links with the largest Jy only. 



For small but finite k, there will be three phases: (i) 
> = for J2/J1 below a critical value {J2/ Ji)ci < 1; (h) 
L = for J2/J1 above another critical value {J2/ Ji)c2 > 
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1; and (iii) coexisting Qi 7^ and (52 7^ for intermediate 
values of J2I J\- 

Longer loops, flux expulsion. The terms of order 
{PP'^)^ represent loops containing up to 4 links. For a 
single (e.g., nearest- neighbor) nonzero exchange constant 
J, PQ^ = 2PPt/J and the energy, to this order, is 



MN 



Tr 



(PPt)' 



Trl 



(10) 



For a fixed "norm" Tr(PPt) = 



[Tr(PPt)]^ 

lower en- 
ergy means larger Tr(PPtppl^). This can be achieved 
by tuning the phases of link fields Qij , as the contribu- 
tions to Tr(PP'''PPt) from loop abcda and its reverse 
are PabP^.PcdPL + Cc - 2Scos$, cf. Eq. ©. Clearly, 
for fixed magnitudes \Pij\, the trace is maximized — and 
energy is minimized — when the flux $ vanishes. 

This establishes the principle of flux expulsion for the 
shortest nontrivial loops (length 4). As the loop expan- 
sion at small k, is organised by loop length, this prin- 
ciple provides the correct ground state as the (ideally, 
only) one in which all loops up to a certain length con- 
tain no flux. For example, it uniquely determines the 
Sp(iV) ground states observed for the uniform triangular 
and kagome antiferromagnetsji 

More generally, we can formulate a conjecture on the 
behavior of longer loops. It provides all gauge-invariant 
information about the phases of link variables Qij . 

Conjecture (flux expulsion): In the ground 
state, the flux $ is zero through all closed 
loops of even length, provided such a fluxless 
state is possible. 

Keep in mind that the tendency to expel flux does not 
always guarantee the actual absence of flux. If the lattice 
is not bipartite, fluxes may be frustrated and will not be 
expelled from every loop. This happens already for the 
triangular and kagome cases and for the latter this makes 
the selection more delicate than believed previously.-^ It 
happens with a vengeance on the pyrochlore lattice. 

Pyrochlore I: the single tetrahedron. The short- 
est loop (of even length) on the pyrochlore lattice con- 
tains four links and is confined to a single tetrahedron. 
There are, in fact, three such loops on every tetrahedron 
and it can be verified that the sum of their fluxes equals tt 
(unless some link amplitudes vanish). The fluxes are thus 
frustrated and cannot be expelled from all three loops. 
Unlike in all previously studied systems, the principle of 
flux expulsion does not point to a unique ground state. 

In fact, for the single tetrahedron, we find a two- 
parameter family of ground states, all with exactly the 
same mean- field energy £'/(JiV A/") = —k,{k,-\-\)/2. These 
have link variables 



Q12 = Q34 = V '^('^ + 1) sin0, 

<9i3 = Q24 = \/ k{k + l)(cos0 — i cos 9 sin (j>) , (11) 

(9i4 = Q23 = V k(k + l)(cos6'cos(/) — i sin0). 




FIG. 1: The two parameter family of ground states of a single 
tetrahedron. The four spins combine to give zero total spin. 



This is quite remarkable as one of the charms of Sp{N) 
is its capacity to yield unique quantum disordered states 
at small n — indeed, this is the first counterexample! Re- 
markably, there is a mapping between these ground states 
and those of classical Heisenberg spins on a tetrahe- 
dron, which can be constructed by (a) parametrizing the 
ground states as shown in Fig. ^ (b) representing the 
spins Si by two-component spinors ipia and (c) translat- 
ing the spinors into link variables Qij oc Eapi^iaipjfi- 

Possible orders. All of these mean- field ground 
states violate a symmetry of the Hamiltonian (Q: a 
point-group symmetry, time reversal, or both. These 
symmetry breakings are best illustrated by the sets of 
states which break only a single symmetry. 

First, breaking the symmetry group of the tetrahedron 
Td, are three bond-ordered states with maximally inhomo- 
geneous bond amplitudes, e.g. Q12 — Q34 = and Q13 = 
Q24 = Qi4 = Q23 7^ 0. These we call the "coUinear" 
states, as their classical counterparts have collinear spins 
[Fig. |2i;a)]. The flux through loop 13241 vanishes; the 
other two fluxes are ill-defined. The valence-bond or- 
der parameter characterizing the broken symmetry is de- 
scribed in Refs. [TJ andlT^ 

Second are states which leave the spatial symmetry 
intact but break the time-reversal symmtery. These have 
6* = i arccos (— 1/3), (j) = ±7r/4. They distribute the flux 
TT equally between the three loops, each receiving either 
-f7r/3 or — 7r/3. The classical analogs of the two states 
have spins pointing at equal angles of arccos (—1/3) w 
109° to each other [Fig. |2Ic)]. The order parameter is 
spin chirality X = (Sa' (Sf, x Sc)}, where a6c is an oriented 
face of the tetrahedron. ^^■^■^ 

Pyrochlore II: the full lattice. The pyrochlore lat- 
tice is a network of corner-sharing tetrahedra. In the loop 
expansion up to O(av^), the system behaves as if it were 
made up of disjoint simplices: the energy is minimized 
as long as each tetrahedron is in any of its ground states 
(|ll|l . This extremely large degeneracy will be lifted, at 
least partially, at 0{k^), which includes loops enclosing 
the hexagons of the pyrochlore lattice. In the spirit of 
degenerate perturbation theory, we must minimize the 
terms 0{k^) among all such possible states. This is a 
problem of considerable complexity. 
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(a) 



(b) 




FIG. 2: Trial ground states of the single tetrahedron: (a) 
coUinear, (b) coplanar, (c) anticoUinear. 



On a first approacli we have minimized the energy over 
a restricted set of trial states in which all tetrahedra 
are in the same state Evaluation of the energy 

for a finite cluster shows that, in accordance with our 
earlier conjecture, the energy reaches a minimum when 
the hexagons enclose zero flux; this choice is similar to 
Sachdev's Qi = —Q2 state on the kagome.i 

We obtain the following energies per site for the ground 
states derived from classical spin states shown in Fig. 12 



p 2 I 0.18750k3 

-^^^ 2 2 1 0.20139^3 Fig.l^c) 



Fig. Ha), 
Fig.Hb), (12) 



The energy is lowest in the state with equal link am- 
plitudes that is related to the classical state with equal 
angles between the spins [Fig. |21[c)] . We have evaluated 
the energy difference between the states (a) and (c) ana- 
lytically and found that {Ei^^) - E(^^))/{JNN) = 
Note that the splitting is rather small: for the nominal 
equivalent of spin 1/2 (k = 1), the energies per spin dif- 
fer by about one percent of the exchange constant. The 
order-from-disorder effect is extremely weak. 

Unfortunately, we have been unable to prove that the 
state described above is indeed the state of lowest en- 
ergy, for we have discovered another state with the same 
energy at 0{k^). Without going into details, we note 
that the other state has a larger unit cell and contains 
tetrahedra with "coplanar" spins [Fig.|2Ib)]. 



Outlook. We have learned that finding the ground 
state of the Sp(iV) pyrochlore antiferromagnet is a hard 
problem. For all previously studied systems, our method 
yields unique ground states in agreement with numerical 
minimizations at lowest non-trivial order in k. In con- 
trast, simplices of the pyrochlore have continuously de- 
generate ground states. Apart from flux expulsion, there 
is no simple principle that can guide the search for the 
groimd state in this case. 

At the same time, our study has produced some use- 
ful insights. First and foremost, we find that there is no 
ground state retaining the full symmetry of the Hamil- 
tonian; such a state is already ruled out at the level of 
a single tetrahedron, in a controlled fashion at small k. 
In S\){N), there is no spin liquid on the pyrochlore lattice 
at zero temperature. Furthermore, there is a huge num- 
ber of nearly degenerate saddle points, which are not re- 
lated by a symmetry, with a splitting 0{k^). These small 
splittings suggest that determining the precise nature of 
the symmetry breaking is going to be very hard and, for 
experimental systems, exquisitely sensitive to small ad- 
ditional terms in the Hamiltonian. 

In the large-iV treatment, tiny energy differences be- 
tween saddle points are made infinite as they come with 
a large prefactor N — > 00. In practice, = 1, and 
therefore local tunneling events will probably play a role. 
Whereas individual saddle points do break symmetries, 
tunneling can lead to symmetry restoration. At the end 
of our investigation, we therefore have to declare our- 
selves agnostic as to the eventual fate of the quantum 
pyrochlore magnet at zero temperature. Whatever or- 
der may be present there will likely melt rapidly at finite 
temperature. 
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